Generation of eigenstates using the phase-estimation algorithm 
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The phase estimation algorithm is so named because it allows the estimation of the eigenvalues 
associated with an operator. However it has been proposed that the algorithm can also be used to 
generate eigenstates. Here we extend this proposal for small quantum systems, identifying the con- 
ditions under which the phase estimation algorithm can successfully generate eigenstates. We then 
propose an implementation scheme based on an ion trap quantum computer. This scheme allows us 
to illustrate two simple examples, one in which the algorithm effectively generates eigenstates, and 
one in which it does not. 
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I. INTRODUCTION 

Since the inception of quantum computation 0|, peo- 
ple in the field have endeavored to find tasks which a 
quantum computer could perform more efficiently than 
a classical computer ||, ||, H, ^. For a detailed intro- 
duction into the field of quantum computation and infor- 
mation, see The algorithm which has by far gener- 
ated the most interest is Shor's factoring algorithm [Q, 
as it enables the cracking of the RSA encryption sys- 
tem 1^. Kitaev j|] generalized Shor's algorithm, show- 
ing how a quantum computer can generate an eigenvalue 
of an arbitrary unitary operator (in the limit of a large 
number of qubits, and not necessarily efficiently). Due 
to experimental difficulties, a large scale quantum com- 
puter (if possible), will not be attainable for a number of 
years. However, small-scale quantum computers are al- 
ready available In this paper, we show how a version 
of the phase estimation algorithm can be implemented 
on a particular 'small-scale' quantum computer, the ion 
trap quantum computer. 

Given some unitary operator U and an approximate 
eigenstate; the goal of the phase estimation algorithm ||, 
is to obtain an eigenvalue of U and leave the quan- 
tum system in the corresponding eigenstate |l^ . To 
accomplish this task, we shall need two quantum sys- 
tems which can be coupled together. One, we shall call 
the index system, the other the target system. The in- 
dex system is initially prepared in the state |0). After 
performing the algorithm, the index system will store an 
eigenvalue of the target system operator, U . 

Traditionally both the target and index systems have 
been qubit registers. In this paper the index system will 
remain a register of qubits, however we shall allow the 
target system to be an arbitrary A^-dimensional quan- 
tum system, where N may be equal to infinity. For a 
more generalized discussion of combining continuous and 
discrete quantum computation see [[l3[. 
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In Sec. II we briefly review the phase estimation al- 
gorithm then derive the analytical results which will al- 
low us to characterize the algorithm's performance when 
using only a small number of qubits. In Sec. [II we de- 
rive the Hamiltonians necessary to investigate the num- 
ber and displacement operators in an ion trap, and con- 
trast the algorithms effectiveness with respect to the two 
different operators. 



II. THE PHASE ESTIMATION ALGORITHM 

In what follows, we shall assume that our index sys- 
tem is a register of m qubits. First, we need to be able 
to perform the operation A([/) on our coupled system. 
A([/) is completely described by defining its action on 
the standard basis states of the index system, coupled to 
an arbitrary target system state. 



A(;7)|j)/|^>, 



Ii®U't\3)i\4')t 



(1) 



where Za/ = {0, 1, 2, . . . , M - 1} and M = 2". As in the 
last line of Eq. (|l|) , we shall continue to omit the subscript 
notation when it is clear whether a ket or operator is 
referring to the target or index system. We begin the 
algorithm by initializing our quantum computer into the 
state 



l*o) = |0)|^). 



(2) 



Performing a 7r/2 rotation of each qubit in the index reg- 
ister results in the state 



A/-1 



(3) 



We now perform K(U) on this state giving 



M-l 



|vI/2)=A(C/)|vI/i) = 1 ^ |^-)[/.|v,). (4) 

The final steps in the algorithm are to perform the uni- 
tary quantum Fourier transform [H on the index register 
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and measure this register However, before applying 
this transform we shaU re- write Eq. (Q) . First we replace 
by it's representation as a sum of eigenvectors of U , 



(5) 



where k sums over the dimensionality of the target sys- 
tem. Hence the state [^2) can be written as 



j=0 k 



Ck\(Pk) 



(6) 



We shall write the eigenvalue associated with \(f)k) as e*'^'= . 
Noting that applied to each eigenvector \(f)k) is simply 
'\(j)k) and changing the order of the summations we 



obtain 



A/-1 



1*2) = 



(7) 



Lastly, for clarity we exchange the order of the systems, 
and replace 4>k with i'KLOk/M^ where S [0, M), 



A/-1 

k ^ j=0 



g27r»ja;./Af|^.^^ (8) 



It is now not hard to show that taking the quantum 
Fourier transform of the index register results in the state 



A/-1 



1^3) = E^^i'^'^-)E/('^fc'j')ij), (9) 

k j=0 



where 

f{uJk,j) = I ^-n(.^) 



1 ^sinU^.) ..(.,-^) 



(10) 



t 1 :wfc=j. 

As we will see shortly, it is helpful to note that 

> |sinc(c^fc- (11) 

for all uJk & [0, M) and j G Zm- A plot of \f{uJk,j)\ is 
shown in Fig. |l| where M = 16 and j has been set to 5. 

Finally, measuring the index register will, with high 
probability, yield an approximate eigenvector. To under- 
stand this, let us begin by looking at the most simplified 
case. Suppose for a moment, that we have ujk G for 
all k. then 



fi^k,j) 

Thus Eq. (^) simplifies to 

^Ck\(l>k)\uJk) 



(12) 



(13) 



If we add the assumption that no two values of k give 
the same ujk (i-C. we have no degeneracy p9|) then upon 
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FIG. 1: A plot of I/I as a function of uj^, with M = 16 and 

i = 5. 

measuring the index register, we will obtain \ujk), and 
hence e"^*-', with probability |cfep, and leave the target 
system in the eigenstate \4>k)- 

Removing the assumption of zero degeneracy, measur- 
ing the index register still allows us to obtain some eigen- 
value e^'^^'il^^ ^ however the target system is now left in 
the state 



(14) 



where fc' = {fc : Wfc = j}, and M = X^fc' |cfc'P is a nor- 
malization constant. 

Finally, we shall remove the assumption that the Wfe 
must be elements of Zm- The probability P{j), of mea- 
suring the index register in some basis state \j) is 



|Cfc/(wfc, j)|^- 



fc 

El 



(15) 



Having measured the index register to be in some state 
\j), the target system is left in the state 



E"" 

fc 



Cfc = 



Cfc/(wfc,i) 



(16) 



whereA/' = X;fc|cfc/('^fc,j)P- 

In order to gain some useful information from Eqs. ( |l5| ) 
and (p^, let us assume that our initial target system 
state IVv is an approximate eigenstate of \(j)q) for some q 
such that 



\{<t>M\ 



(17) 



Remembering that ujq will be some real number between 
and M , we define [w^J to be the nearest m-bit integer 
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less than cuq, and [wg] to be the nearest m-bit integer 
greater than ujg, where modulo M has been assumed. 
The probability of measuring the index register in either 
the state \ [iOq\) or is 

k 



UJa 



> |c,/K,Lc^,J)p + |c,/K,ra;,l)P 

> |c,p 2 sinc2(0.5) 

> 0.8p. (18) 

Hence, with probability greater than 0.8p we will obtain 
an approximate eigenvalue associated with \(j)q) which 
differs in phase from the actual eigenvalue by less than 
27r/2™. Thus, if p is reasonably large, we have a high 
probability of finding the best estimate of the eigenvalue. 
However, as we shall see, large p does not imply that we 
will improve on the approximate eigenstate. 

Suppose we measure the index register in the state 
\[ujq]), where [iVq] denotes the closest m-bit integer to ujq. 
(N.B. This will occur with probability greater than OAp, 
as \f{ujq, [iJ^ijDP > 0.4.) The key question that we wish 
to address in this paper is: has our initial approximate 
eigenstate improved? Letting p' = Ic^p, we are effec- 
tively asking what bounds can be placed on p'l For an 
arbitrary U it is obvious that the upper bound of p' = 1 
can be obtained by setting lijj) = \(j}q). We now inves- 
tigate the lower bound by dividing the eigenstates into 
three disjoint sets, 

s - 

g = {g:g^q,\LOg-[uq]\<l}i,iid (19) 
H = {h: >1}. 



We now have 



p\f{ujq, [uJq] 



(20) 



with 



see 



+ 5^ic,n/(c.,,K])i= 



(21) 



Using Eqs. ( p^ and (|T^), it is not hard to show 

0.4 < \J{Uq,[ujq])\^ <1 

o< |/(c^s,K])l' <i 

0< |/(c^,,K])|2 <A 



(22) 



where A = |/(1.5, 0)p. As m increases A tends to {j^Y 
0.045. However for our analysis it is sufficient to note that 



0.045 < A < 0.05 for m > 3. Eq. (H) leads to the lower 
bound 



P' > 



p|/(a;„K])|2 + (l-A)G + A(l-p) 



(23) 



where G — |cgp. Fig. ^contains a plot of this lower 
bound as a function of G for |/(wg,[a;g])p — 0.6 and 
various values of p. The circles indicate the points at 
which the minimum of p' equals p. Thus, we see that by 
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FIG. 2: The lower bound on p' as a function of G for 
\f(ujq, [oJq\)'\^ = 0.6 and various values of p. The circles in- 
dicate the points at which the minimum of p equals p. 



endeavoring to make G as small as possible, we increase 
the amplitude of |0q). For a given U and G can be 
made arbitrarily small by increasing m. However, we are 
interested in the performance of the algorithm for small 
values of m. We shall now look at G's dependence on 
U and \'ip) by attempting to create eigenvectors for both 
the number and displacement operator in a ion trap. 



III. AN lONTRAP IMPLEMENTATION 



We first derive the Hamiltonian for A(?7), where U is 
the evolution operator associated with the number op- 
erator, and investigate the phase estimation algorithm's 
performance for various initial states. We then derive the 
Hamiltonian for the more complicated case of U being 
the displacement operator. For both of these examples 
the index register will be two electronic levels of m ions 
in a linear ion trap, and the target system will be the 
center-of-mass (CM) vibrational mode of the ions. 
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A. The Number Operator 

Consider the standard Hamiltonian of the one dimen- 
sional harmonic oscillator, 



H = hu){a) a + 



1, 



(24) 



where and a are the creation and annihilation opera- 
tors. Ignoring the over-all phase contribution of the zero 
energy state, the unitary operator we will first be ana- 
lyzing is 



U{t) = e"^" 
In this case, A([/) is given by 



A(f/) 



exp 




where 



(25) 



(26) 



(27) 



The inversion operator for each ion is defined by CTz — 
(|0)(0| - |l>(l|)/2. This Hamiltonian can be obtained for 
interaction times greater than the period of the CM vi- 
brational mode by applying a set of far-detuned standing 



wave pulses to the ion 1 15 



We begin our analysis by initializing the CM mode in 
some phonon number state \n) [ p^ , and setting ujt — 
27r(l — 1/M). It is important to note that we are as- 
suming that all the higher vibrational modes are in the 
vacuum state. Assuming no errors, applying the phase 
estimation algorithm results in the index register being 
measured in the state \n mod M) and the target system 
is left unchanged. If we now let cot be some arbitrary 
value, applying the algorithm will leave the target sys- 
tem unchanged, and the index system will be measured 
in the state \j) with probability 



P{j) 



. LotnM 



(28) 



Let us consider the more interesting situation where 
the target system is initialized in some coherent state 
I a). We can utilize the phase estimation algorithm to 
transform the state of the target system into an approx- 
imation to a Fock state. 

For example, suppose we use four index qubits, tot — 1 
and we choose to approximate the Fock state |n = 9) by 
using the coherent state ja = 3). In this example we per- 
haps might think that ja = 3) is not a good approximate 
state because p « 0.13 however the fact that G < 0.035 
indicates the algorithm should work well. Applying the 
algorithm and measuring the index register in the state 
|9), we obtain p' w 0.93. The initial and final target state 
for this scenario is shown in Fig. |^. 
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FIG. 3: Fock state distributions for the target system initially 
in a coherent state with a = 3, and the state of the system 
after applying the phase estimation algorithm and measuring 
the four index qubits. 



Having shown that the phase estimation algorithm can 
be used to generate Fock states from coherent states, we 
now attempt to generate eigenstates of the displacement 
operator. 



B. The Displacement Operator 

The displacement operator applied to the CM vibra- 
tional mode is defined as 

D{a) = exp(aa'f - a*a). (29) 

Thus the operator we wish to apply is 

m — 1 



AiD) = exp I ^ iJ J , 



(30) 



3=0 



where the Hi are now defined as 



Hj = ih[aa) - a*a)2J'(CTp) + i). (31) 

It has already been shown ||l^ that conditional displace- 
ment operations such as the Hamiltonian in Eq. (|3^) can 
be performed in an ion trap. 
It is not hard to show that 

D(de'(^+^))|a,e) w e^^^l^l^^^a, e) (32) 

for large values of squeezing parameter r and where a = 
\a\e"t',e^ re^^'t' and 



\a,e) = S{e)Dia)\0) 



(33) 



is a squeezed coherent state. Thus the squeezed coherent 
states \a,e) form approximate eigenvectors of the dis- 
placement operator Z3(de*''^''"^)). 
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Without loss of generality we can set cj) = in which 
case the eigenstates of the displacement operator are 
simply the position eigenstates. It is then not hard to 
show that for small fixed m, G « 1 — p, which leads to 
p' w p. Thus applying the phase estimation algorithm 
to squeezed displaced states does not produce improved 
eigenstates of the displacement operator. 



implementation could be performed using cavity QED, 
allowing generation of photon number states with only 
small numbers of trapped atoms. We have also shown 
that the algorithm's performance depends on the relation 
between the approximate eigenstate and the spectrum of 
the operator. We can gauge the algorithm's performance 
by calculating a parameter G. 



IV. CONCLUSION 

We have shown that the phase estimation algorithm 
can be used to generate eigenstates of the number oper- 
ator, even when we severely limit the size of the index 
system. It would be interesting to see if an analogous 
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